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Abstract 

In a large variety of spectroscopical applications Bloch-Boltzmann equations (BBE) play an 
essential role. They describe the evolution of the reduced density operator of an active atom 
which is coupled to radiation (Bloch part) and which interacts collisionally with the perturber gas 
(Boltzmann part). The standard approach to the collisional part is well-known from the literature. 
It preserves hermiticity and normalization, but the question whether it preserves positivity seems to 
remain open. The completely positive BBE were recently derived via the general master equation 
techniques. These two approaches are applied for a model of n-level nondegenerate atom. We show 
that within this model both approaches to the collisional part of BBE are equivalent - give the 
same physical predictions. The approach based upon master equation techniques guarantees the 
preservation of hermiticity, normalization and positivity. The proven equivalence ascertains that 
the standard approach also preserves positivity. Moreover, some aspects of the standard derivation 
(which atomic states do contribute to the evolution) are clarified by the established equivalence. 
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I. INTRODUCTION 



One of typical, spectroscopically important experimental situation occurs when a gaseous 
mixture is irradiated by the external light source. The mixture consists of active atoms which 
couple to the incident radiation and of a usually much denser inert gas. The perturbers atoms 
collide with the active ones, thereby influencing their behavior. Then, many properties of 
such system can be studied, both experimentally and theoretically. The literature devoted 
to such problems is large, so we indicate only some essential monographs 

Usually, only the active atoms are of interest, hence they must be theoretically described 
within the density operator formalism. The equations of motion for the active-atom density 
operator may be called, in the absence of any better name, Bloch-Boltzmann equations 
(BBE). The Bloch part describes the interaction between the active atoms and radiation. 
It is a generalization of the well-known two-level optical Bloch equations to a more general 
multilevel case. The coUisional interaction between active atoms and perturbers is given by 
suitably constructed collision integrals. This contribution to equations of motion might be 
called the Boltzmann part. 

The recent advances in the fundamentals of quantum mechanics and in the quantum 
information theory have shown the importance of the preservation of the basic properties of 
any density operator: hermiticity, normalization and positivity (for an excellent review, see 

and the references given therein). These essential properties of the density operator must 
be preserved by any theoretical formalism. The aim of this work is to discuss this point 
for the spectroscopically important situation which was sketched above. Since radiative (or 
Bloch) part of the corresponding equations of motion is already well investigated we shall 
focus our attention on the coUisional (Boltzmann) part. 

n 

In the recent paper fd\ a gaseous mixture of two species: A - active atoms and P - 
perturbers, was considered. The densities of these two components are assumed to satisfy 
the relation: Na -C Np, and an equation of motion for the density operator of the A-atom 
interacting collisionally with the perturbers, i.e., is rederived. This approach is based upon 
general master equation (ME) techniques (for a review, see ^|) in the spirit of Lindblad- 
Gorini-Kossakowski-Sudarshan method. The employed technique ensures that the A-atom 
density operator possesses all the necessary properties: it is hermitian, normalized and 
positive-definite for all instants of time. It is perhaps worth mentioning that the Bloch part 
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(i.e., radiative one) of the equations of motion for the active-atom density operator is usually 
derived within ME techniques (see Therefore, this contribution to BBE is certain to 
preserve the mentioned properties of the density operator. This is also the reason why we 
restrict our attention to the collisional part of BBE. 

The other approach to the derivation of the Boltzmann part of the BBE is known since 



the pioneering work of Snider 



. Then, it was refined by other authors and employed in a 



variety of practical applications (see, for example 0, [l^). The monograph P] by Rautian 
and Shalagin seems to give the most comprehensive review of the derivation of collision 
integrals appearing in BBE, especially for specific, spectroscopy-oriented applications. In 
the following, we shall call their presentation the standard approach to the derivation of 
BBE. It is, perhaps, worth noting that the derivation given by Rautian and Shalagin leads 
to the appearance of some Kronecker-like delta factors which perform the role of "state 
selectors" - select the atomic states which contribute to the evolution of the A-atom density 

n 

operator. The authors themselves say that their procedure is open to question (.2], p. 42). 

Working within the standard approach it is relatively straightforward to prove that it 
preserves hermiticity of the atomic density operator. The proof that the normalization 
is also preserved requires one to invoke the optical theorem of the quantum multichannel 
scattering theory. However, we are not aware of any studies in which the preservation 
of positivity is investigated. Discussion of this point is the main aim of our work. The 
secondary aim of our work is to investigate the validity od the " state selection" mechanism 
proposed by Rautian and Shalagin 

It is not our aim to present the details of the derivations, or the underlying physical 
and mathematical assumptions, of the two indicated approaches to Boltzmann part of BBE. 
These aspects are well-documented in the literature 3,1^, so it seems to be no need to repeat 
them here. We shall only use the results relevant to the main subject of our discussion. 

Sec. II is devoted to brief presentation of the collisional contributions to BBE obtained 
within the master equation approach and in the standard one. We do not derive them, but 
simply state the results - the results which are essential for further discussion. 

In the next section we adapt the general formulas of the previous section to the model 
of n-level atom with nondegenerate energies. We construct the collisional parts of BBE 
corresponding to both approaches. In the last subsection we argue that the obtained results 
are, in fact, equivalent. This allows us to state that the standard approach (within the 
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adopted model) preserves the positivity of the active-atom density operator. The proof of 
this fact constitutes the main result of this work. 

Finally, the fourth section contains some remarks which may be useful for some further 
research, that is for discussing the preservation of positivity for more general models studied 
within the standard approach (master equation approach is guaranteed to do so). Moreover, 
we hope that some of our remarks will be useful to provide the standard approach with 
some more rigorous standing. Namely, the equivalence of both approaches validates the 
the Kronecker delta- like "state selective" factors which appear in the standard approach of 
Rautian and Shalagin. 



II. TWO APPROACHES TO DERIVATION OF THE COLLISIONAL PART OF 
BBE 

A. Master equation approach 

The derivation of the Boltzmann part of BBE via the master equation technique is given 
in the recent paper jof, where the necessary assumptions are also discussed. This is a 
mathematically rigorous, although fairly involved method. It is not our purpose to obscure 
the physical discussion by mathematical technicalities, therefore, we will present here only 
the most essential results of the theory given in Qj. 

In order to ensure the preservation of hermiticity, normalization and positivity of the 
reduced atomic density operator, the corresponding equation of motion must be of the 
following general (so-called Lindblad-Kossakowski-Gorini) form (see also \^) 

- ^iBapa + PaBa), (l) 

where pa is a (reduced) density operator of an A-atom. It is parameterized by an index 



a which, under the suitable additional assumptions 



3, 



can be shown to correspond to the 



velocity of a considered atom. Let Tia denote a Hilbert space of the atomic states. Then 
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the quantities appearing in Eq. (jT)) are defined as mappings (operators): 

Ha = H}^ '■ — * T^a, (Hamiltonian) , (2a) 
S^^ : Hp ^ Ha, (2b) 
{S^J-.n^^Hp, (H.C.of^i,), (2c) 

^« = EE(^i/^l' (2d) 

The general master equation ((T)) is then adapted to our needs - to describe the active-atom- 
perturber coUisional interaction. These steps are also discussed in j^l, so we only state the 
results essential to further discussion. 

We consider a ra-level A-atom immersed in much denser perturbers. The perturbers 
(assumed to be structureless particles) thermalize very rapidly, hence their distribution is 
simply Maxwellian 

W^^^)(v) = fAV'expf-^V (3) 



with Up = 2kBT/mp, where rUp is the mass of the perturber atom. 
Let us now take the Hamiltonian of the free A-atom as 

n 

HA = Y,nuJk \k){k\ (4) 

k 

where the eigenfrequencies ujk may, in general, be degenerate. 

Next, let {Sa\ be a basis in the space of operators acting on the Hilbert space of A-atom 
states (I A;)}. These operators satisfy the relation 

[Ha, Sa] = hQaSa, a =1,2, , (5) 

where the quantities Qa are identified as Bohr frequencies. 



Within this framework, the coUisional part of the master equation becomes 



d \ 1 r 

P(v) ,= -7;^^ba{^) SlSb, p(v) 



dt \coii. 2 

a,b 



(+) 



+ J2jdY' /C,fe(v ^ w')SaPi^')Sl, (6) 

a,b 

where the (-I-) subscript denotes the anticommutator, and p(v) = p(r, v, t) is the reduced 
density operator of an A-atom with respect to internal variables (states) but a phase-space 
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distribution with respect to position and velocity. The relaxation (collisional) rate 7;,a(v) is 
specified as 

76a = 76a(v) = / rfv ' /Cba(v' ^ v). (7) 



Finally, it can be shown 0, Q] that the matrix lCab{^ ^ v') is expressed as 

= 2iVp 5^^,n, j d^ri j d^r H^(^Hv' - v,i) 



X 5' - ^ 



V - V IV,, - V,.i 



a 



X 6 iv"^ -v"^, + 



2hn, 

rl "I 

X /a(Vr <- V^rl) /fe (^r ^ V,,i). (8) 

The employed notation is as follows, rria is the mass of an A-atom, while fi is the reduced 
mass of A-P colliding partners, v^ and v^i are the relative velocities. Np is the density 
of the perturber gas. The functions /a(v,. ^ v,.i) are connected with usual (taken in the 
center-of-mass frame) scattering amplitudes: 

n 

J2 - Vrl)^a = /(j, V, ^ k, V,,i) \j) {k\ . (9) 

a j,k=l 

We note that the K,ab{^ ^ v') matrix is clearly hermitian and positive definite. Hermiticity 
of the matrix Kabiy ^ v') implies that 7^a(v) = 7afe(v). These ensures the preservation of 
hermiticity of the A-atom reduced density operator. 

The factor 5q,^^q,^ in Eq. (jH)) has the sense of the Kronecker delta 

Jo for n^^n,, 

(>na,n, = < (10) 
I 1 for Qa = 

This factor appears due to the secular approximation which is necessary in derivation of 
master equation (as is clearly shown in |^]). The given /Caf,(v ^ v') matrix also ensures the 
momentum and energy conservation. 

We stress that the resulting collisional ME (jHl) preserves all the necessary properties of the 
A-atom reduced density operator p(v). Preservation of hermiticity follows from hermiticity 
of /Cafe(v ^ v') matrix. Next, relation ((7j) ensures that 

d 



dt 



dw Tr{p(v)} = 0, (11) 



as necessary for preservation of normalization. Finally, preservation of positivity is ensured 
by general Lindblad structure of the master equation (jH)). 



B. Standard approach 

The standard derivation of the Boltzmann part of BBE as presented by Rautian and 
Shalagin [2] is rather lengthy and fairly complicated. It is based upon two physical assump- 
tions: Na <^ Np, so that the A-P collisions are binary, and the duration of the collision is 
by far the shortest time scale, so that the impact approximation is valid (see j^, p. 31). The 
general von Neumann equation for the density operator for the entire system is truncated 
(traced) to an equation for a single A-atom. The interaction with the perturbers is then 
considered within the framework of the time-dependent scattering theory. The collision in- 
tegrals are then expressed (similarly as in [8j) in terms of the elements of the scattering 
T-matrix, which are subsequently reexpressed by usual scattering amplitudes. Further steps 
consist in semiclassical approximation which leads to the Boltzmann terms of the following 
shape 

d . 

Paa'(v) „ = 
coll. 



dt 



r(aa',v I aiai')Paia;(v) 
+ y rfvi /C(aa', V I aia'i,vi)p^^^/(vi), (12) 

where p^^ > (v) denotes the matrix elements of the A-atom reduced density operator (which 
has the same sense as in the ME approach). The indices a should be understood as multiple 
ones (atomic states may be labelled by several quantum numbers). 
The collision rate appearing in the first term is given as 



Y{aa , V I aio;^) 



y rfv, iy(^)(v- V,) 



- /*(«', V,. ^ a'l, v^)5(Wc.'a;)'^a«i • (13) 

We note that this collisional rate is given by the elastic forward scattering amplitudes. S^ai 
are simple Kronecker-type deltas, while the factors 6{LJaai) = ^{^a — ^ai) have meaning 



similar to that defined in Eq. (jlOj) . These delta factors ensure energy conservation. Their 
origin and significance will be discussed later. 

The second term of Eq. (|12j) contains the collision kernel specified as 



/C(aa ,v I aia^jVi) = 

X 2Np j dv,, j dvri W^^\vi - v^) 



X 5 ( V - Vi - — (Vr - ^rl, 

rrir 



X /(a,Vr ^ ttl, Vrl)/*(tt', <- V,.i) 

X 6 (^v^, - v',, + ^{E^ - E,J^ . (14) 

The peculiar feature of the Rautian and Shalagin derivation consists in the appearance 
of the Kronecker-type delta factors. These factors play a selective role, indicating that not 
all matrix elements p^.^^,' (in the rhs of Eq. ()12|)) contribute to the evolution of p^a'- The 
origin of these factors is explained by Rautian and Shalagin in the following way. 

The T-matrix elements ( a, | T | ai, p^i ) (taken in the center-of-mass frame) include 



phase factors of the type of 

'it 

exp 



n 



(15) 

When Bohr frequencies ujaai = {Ea — Ea-^)/h are non-zero, the corresponding exponentials 



oscillate rapidly and their contribution to the overall evolution averages out virtually to zero. 
In other words, only those states for which Bohr frequencies are close to zero contribute 
significantly to collision integrals. This argument of Rautian and Shalagin gives rise to the 
5(^^aoi)-type factors in the collision rate and kernel. It must be, however, stressed that these 
simple delta-hke terms appear due to the assumption that the perturbers are unpolarized. 
If this assumption does not hold, the structure of the corresponding delta-like factors would 
be different and more complicated since these factors would include Bohr frequencies also 
for perturbers. 

Rautian and Shalagin discuss the role of the exponential factors but do not carry their 
calculation as far as we did They retain factors (jlSp in their formulas and comment 
only verbally on their significance. Moreover, up to our knowledge, no other authors conduct 
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such a discussion. The reason, perhaps, is that even Rautian and Shalagin doubt upon the 
vahdity of such a "selective mechanism" expressed by the exponential phase factors 
and consequently by delta-type factors (see their discussion at the top of p. 42 in 
The significance of the discussed delta-like factors reduces to, roughly speaking, that "like 
is excited by like" p. 42). This means that the equation of motion (jl2|) connects 
populations with populations and coherences with coherences. 

In the forthcoming we shall return to the discussion of this point, when we shall compare 
the ME results with the standard ones for a more specific model of an active atom. This will 
allow us to shed some new light onto the role played by the "selective" delta-like factors. 

Furthermore, we note that it is relatively easy to show that 

r*(aa , V I aia^) = r(^a a, v | a^ai), (16) 

and similarly 

}C*{aa , V I aia^, vi) = IC{a a, v | a^^ai, vi). (17) 

Both these relations ensure that the evolution given by Eq. (|T^ preserves the hermiticity 
of the A-atom density operator. 

The preservation of the proper normalization of the density operator on one hand follows 
directly from the general formalism employed by Rautian and Shalagin. On the other hand, 
it should be also possible to prove that the summation and integration over v of the diagonal 
equations (|T^ yields zero, as required by normalization. The complicated structures of the 
rate F and kernel /C make it a rather formidable task (at least in general). It seems that it 
is better to use the general equations for some specific model of A-atom structure. Then, 
checking that normalization is indeed preserved should be much easier. We shall do so in 
further section. It is, however, not clear whether the positivity of p(r, v, t) is also preserved. 
There seems to be no compelling, rigorous argument to state so for certain. We shall later 
return to the discussion of this very important issue. 

Finally, we note that the collisional rate (fT^ has the following interesting property 



Re 



2 

which will be useful in the further discussion. 



F(aa , V I aa ) 

^r(^aa, V I aa) + r(^aa, v [ a a ) , (18) 
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III. DISCUSSION 



A. General comments and atomic model 

Equations (jH)) and (jl2j) representing two approaches to the QMBE are of the similar, 
although not necessarily identical form. These external dissimilarities have led us (see 1^) to 
the supposition that the standard approach may not preserve the positive definiteness of the 
A-atom density operator. This supposition was somewhat strengthened by two additional 
facts. 

The general structure of the collision kernels: (jH} In the ME approach and (jl4|) in the 
standard one, is quite the same. The only external difference consists in the structure of 
the "state selective" delta- like factors. At the first sight it is not at all clear whether these 
factors lead to the same "state selection" mechanisms. 

Secondly, the collision rate F is given in Eq. (fT^ by the difference of forward scattering 
amplitudes, while ■jba defined in clearly contains products of scattering amplitudes, as it 
follows after insertion of (jH)) into ((Tj). 

These arguments seem to support the supposition that the standard derivation is not 
certain to preserve the positivity of the A-atom density operator. To clarify these points 
we shall consider the A-aXom with multilevel but nondegenerate structure. Hence, we once 
again take the free atom hamiltonian as: = XIa; \k) {k\, with kets | k ) constituting 
an orthonormal and complete basis in the Hilbert space of atomic states. The free evolution 
of the elements of the atomic density operator is given as 



while we assume that Bohr frequencies Umn 7^ ^jk for different pairs of indices. 
B. Master equation approach 

The general structure of the master equation is given in Eq. (jU)) and it must now be 
adapted to the presently considered model of multilevel nondegenerate atom. The choice of 
the operator basis is in this case obvious. We simply take 




(19) 



Sa < ^ Pjk = \j) {k 



(20) 
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Hence, index a used previously to enumerate the operator basis is now replaced by a pair 
of numbers {j,k). Moreover, Bohr frequency Qa corresponds now to Ujk- The considered 
density operator can then be expanded in the chosen basis as 

P(v) = E/^^'^(^)^^-^- (21) 

Using the adopted identifications all the terms in Eq.® can easily be computed. Since we 
are mainly interested in the comparison of the two variants of the collisional terms of BBE 
we shall omit the computational technicalities. The kernels JCab = ^jk,mn (and consequently 
the rates = 'yjk,mn) contain the "state selective" factors Sq^^q^ = 5{ujjk — ojmn) acting as 
Kronecker deltas. Careful but simple computation of all the necessary sums leads to the 
following collisional equations of motion: for populations we get 

at \coii. 

+ ^ 1 dv' K,mk,mk{^)pkk{^'), (22) 

k 

where we notice that the presence of the summation reflects the fact that inelastic collisions 
are also accounted for. On the other hand, for coherences we get 



Pmn(v) 



coll. 



= - ^[7mm(v) +7rm(v)]pm„(v') 

mm,nn 

(v ^ v')pm„(v ). (23) 
In the two just given equations we have introduced a convenient abbreviation 

7mm ( Vl ) = ^ 7im,jm ( Vi ) 
j 

= ^ y JCjmJmi^ ^ Vi). (24) 
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The collision kernel appearing here is now written in the following form 



^mj,nk{^ ^ Vi) — 

= 2Np j j dwr, iy(^)(vi - V,,: 
X 5^ 



V — Vi 



2 J^ll) jYi, j 



(25) 



Once again we feel it necessary to stress that the resulting equations (j^ - ip!^ preserve 
hermiticity, normalization and positivity of the atomic density operator p(v). 



C. Application of the standard approach 

Equation ()12|) has now to be transformed to suit the currently investigated model. Gen- 
eral (multi) indices a, a ' should be replaced by numbers j, k, etc. The corresponding changes 
are then to be made in the expressions (fT!?|) and ()14j) where the summations involving the 
"state selective" delta-like factors can now be easily performed. This leads to the following 
equation of motion for populations 

d ^ \ 

<jt \coll. 

+ ^ dviJ[mm,\\kk,vi)pkk{^i), (26) 
k 

where we again see the contributions from inelastic collisions. The corresponding equation 
for coherences reads 



^ Pmn(vJ 



(m^n) 
coll. 



+ j d\i J[mn,\\mn,\i)pmn{^i)- (27) 
The collision rate ()13|) transformed to suit the presently considered model is now given as 



^mniy) =T[mn,w\mn) 
'2'Kh 



Nf 



^ j dWr Vr(^)(v- v.) 



I p 

X [/(m, Wr ^ m, \r) - f*{n, <— n, v^)] . (28) 
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and the collision kernel is of the following form 

j{mn,Y\jk,\i) =2Np j dwr j dw^i W^^\wi-Wri) 



X (5 ( V - Vi - — (Vr - Vrl) 

nia 
2 



X /(m, ^ j, Vri) f*{n, V,. ^ k, v^). (29) 

Definition implies that rmn(v) = r^ml"^)- Moreover, general relation (fT7|) yields 
JT* (mn, v|jfc, vi) = J7'(nm, v|A;j, vi). These two facts ensure that Eqs. and 
preserve hermiticity of the density operator. 

To prove that the normalization is retained properly, we need to show that relation (jlip 
is satisfied. Due to Eq. ^I^i is equivalent to the condition 

rfcfc(vi) = ^ / dv j{mm,v\kk,vi). (30) 

m 

Obviously, definition of the collision rate implies that 

Tkkiy) = Np (^^^ I d^r W^(^)(v - V.) 

xIm{/(A;,v, ^fc,v,)}- (31) 
Since the discussed formalism allows for inelastic scattering, we need to use multichannel 
scattering theory 12] . Optical theorem allows us to express the imaginary part of the elastic 
forward scattering amplitude by the total cross section o"r(fc, v^) for scattering from the state 
\k) ,'Vr. Thus, we cast the Ihs of (jHUj) into the form 

rfc,.(vi) = Npj dvr l^(^)(vi - V,) |v,| arik, v,). (32) 
The total cross section can be written as a sum 

MA:,v.)=X:/^^(--) (33) 

where in the rhs we have differential cross sections corresponding to scattering from state 
\k) ,'Vr to \m) , Vj-i and where the integration is performed over the angles specified by the 
direction of final velocity. Then, the collision rate ()H2|1 takes the form 

rfcfc(vi) =Np'^ J d\r j dn{\r 



'rl. 



l^(^)(v,-v.)|v.|-^^. (34) 
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On the other hand, rhs of Eq. contains square moduh of scattering amphtudes (as it 
follows from ((221) )• From multichannel scattering theory we have 

\f{m,Vr ^ k,Vri)\ = —rj . • (35) 

(note the reversed roles of relative velocities and Vri). Inserting the kernel (j^^ into the 
rhs of (j30|) and using (j35|) we can perform all the necessary integrations. Then, we arrive 
at the expression identical with rhs of (jH^ . This completes the proof of relation (jHIIl) and, 
therefore, we see that the considered model ensures preservation of the proper normalization 
of y4-atom density operator. We note, however, that the question of positivity preservation 
still remains open. 

The general property ()18|) of the collisional rate in standard approach allows us to write 
for the presently considered case 

Re { r^„ ( v) } = ^ ( ^mm ( v) + r„„ ( v) ) , (36) 

which has important consequences. Eq. ()27|) describes the collisional evolution of coherences. 
The term containing imaginary part of Tmn can be written separately. Then it can be 
combined with hamiltonian (unitary) part of the evolution and identified as the coUisionally 
induced atomic frequency shift. Therefore, only the term containing RelFmn} contributes 
to relaxation part of Eq. ()27|). which is therefore replaced by the following equation 



at 



coll. 



(r^m(v) + r„„(v)jp„„(v) 

+ j dvi J{mn,v\mn,vi)pmni^i), (37) 



Moreover, we note that within the discussed model only the rates Tkk{'^) given in (jnn|) - (jHT|) 
or in ()34|) are of importance. 

Summarizing, we can say that within the standard approach the evolution of the density 
operator of a multilevel nondegenerate atom is governed by Eqs. (pUj) and (jHTj) for populations 
and coherences, respectively. Finally, we note that Rautian and Shalagin give correspond- 
ing equations which do not fully agree with the above given results |l3i. We attribute this 
discrepancy most probably to misprints. Moreover, we have fully used the "state selec- 
tive" delta-like factors which, unfortunately, are not consequently resolved by Rautian and 
Shalagin who retain exponential factors like one in p5|) . 
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D. Equivalence of both approaches 



Collisional equations of motion for both discussed approaches are given by formulas 
and for populations, and by Eqs. (f2^ and (jHTj) for coherences, respectively. By in- 
spection, we see that these two pairs of equations have the same formal structure. On the 
other hand, comparing Eqs. ()25p and ()29p giving the kernels for respective approaches, we 
see that they are, in fact, identical 

j{mn, v\jk, vi) = /C.mj>fc(v ^ vi). (38) 

Next, from Eq. (gH) it follows that 

am(vi) = ^ j d\i ]Cjm,jm{^i ^ v) 



7n 



^ / rfvi j{jj,Yi\mm,Y) 



= rfcfc(vi), (39) 

where, in the last step, we have used relation ()3Up which was proved for the standard 
approach. 

We conclude, that the pairs of equations and pHj) - (p?7jl are not only of the 

same formal structure but are strictly identical. This allows us to answer the question on 
the preservation of positivity of the atomic density operator within the standard approach. 
Since master equation technique is guaranteed to do so and the standard approach yields 
the same results, it also preserves all the necessary properties of the atomic density opera- 
tor. Obviously, this conclusion is valid for the considered model - a multilevel atom with 
nondegenerate levels satisfying the requirement that Ujk 7^ ujmn for two different pairs of 
indices. 

It is perhaps worth noting that proving the equivalence of both approaches we have used 
Eq. ()37|) for evolution of coherences in the standard approach. Certainly, this equation is 
valid due to the possibility of replacing r^n in Eq- by the right hand side of Eq. (jHBjl . 
Hence, the latter one is essential in the proof of equivalence and it follows directly from 

n 

definition (j2Hl), see also ^13.]). 
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IV. FINAL REMARKS 



We have compared two different approaches to the derivation of the colhsional (Boltz- 
mann) part of the spectroscopically important Bloch-Boltzmann equations: the master equa- 
tion approach and the standard approach as reviewed by Rautian and Shalagin P]. We 
have shown that both approaches for a multilevel atom with nondegenerate levels and with 
nondegenerate Bohr frequencies are equivalent. Within the given model both preserve the 
fundamental properties od the atomic density operator: hermiticity, normalization and pos- 
itivity. We have thus proved that the supposition stated at the beginning of Sec lIIII was 
wrong and both approaches, at least within the discussed model, are equivalent. 

As it is seen from our comments lifting any of our assumptions may lead to different 
conclusions. First of all, we note that allowing for overlapping line profiles, that is for ujk 
not necessarily different from (with {j,k) ^ {m,n)) may lead to other results. In the 
above considered case the colhsional equations of motion for coherences (PHj) and similarly 
dnZj) connect the given coherence pmn only to itself. This is clearly due to the assumption 
that Ujk 7^ t^mn- In the above presented model the populations (see Eqs. or 
are coupled to other populations. This population transfer is obviously due to inelastic 
collisions which induce excitation-deexcitation processes. It may be expected that when 
the assumption ujjk 7^ oOmn is not valid, the same would happen to coherences - inelastic 
collisions would induce polarization transfer, that is would couple coherence pmn to other 
coherences. However, details of such couplings would certainly depend on the degeneracies 
between different Bohr frequencies. It is, therefore, difficult to give any predictions on 
the specific couplings. Perhaps, the best way is to examine some concrete model with the 
discussed degeneracies and thus exhibiting polarization transfer. Then, the problem whether 
the standard approach to BBE still retains all the required properties of the atomic density 
operator needs to be reexamined. In this context we may say that the master equation 
technique is advantageous. Due to its rigorous mathematical background it will certainly 
preserve all the necessary properties of p(r, v, t) . Indicated problems are clearly of interest 
and seem to be a good subject for further investigations. The question of conservation of 
positivity could be answered again and the advantages of both approaches would be weighted 
again. 

Secondly, it is important to remember that we have taken the perturbers to be structure- 
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less. It is, however, not an essential simplification. It is straightforward to generalize our 
results to unpolarized perturbers, as it is consequently done by Rautian and Shalagin. Ad- 
ditional degree of freedom (that is /? - an index indicating an internal state of the perturber) 
will result in additional summations in the expressions for the collision kernels and in addi- 
tional terms in the Dirac deltas responsible for energy conservation. The "state selective" 
factors either in the standard or in the master equation method would remain unchanged. 
Moreover, it seems that the assumptions concerning the perturbers are less constraining. 
The reason, having purely physical but not mathematical background, seems to be simple. 
In the majority of spectroscopically interesting situations the noble gas serves as perturbers. 
The excitation energy of noble gas atoms is usually beyond the region of interesting ener- 
gies of A-atom transitions. So, the perturbers indeed act as structureless particles. As we 
already mentioned, it does not seem difficult to generalize our master equation approach to 
perturbers with full internal structure (discarding even the assumption about nonpolariz- 
ability). The question is whether it is experimentally interesting or relevant. But on the 
other hand, the problem of positivity conservation within the standard approach would need 
to be reexamined again. 

Rautian and Shalagin introduce the "state selection" via the exponential phase factors 
()15|). They also express their reservations about the validity of such an approach. It seems 
that, at least within the model discussed in this work, their method is equivalent to the 
secular approximation inherent in the master equation. The proved equivalence of two 
approaches seems to clarify the relationship between two derivations of collisional terms 
in the BBE and to strengthen the arguments given by Rautian and Shalagin. However, 
when perturbers are allowed to be polarizable the phase factors of Rautian and Shalagin 
do not necessarily coincide with the delta factors resulting from secular approximation in 
the master equation approach. Consequences of this, also present another problem which 
deserves further investigations. 

We hope that our work is useful to clarify some questions concerning the various ap- 
proaches to the derivation of Boltzmann parts of BBE. We also hope that the results given 
here would be a useful starting point for some further research. 
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p. 46, Eq. (2.157)) is written as follows 



r(i^)(v) = Np (f^) I dv. Ty(^)(v - V.) 

X [f{m,Vr ^ n,Wr) - f*{m,Wr ^ n,Wr)\. 

X Im[/(m, ^ n,Wr)\, 

with notation changed to one used in this work. We note that the rate Fmn given by Rautian 
and Shalagin is, for populations, the same as our one given in Eq. (|31j) . However, for coherences 
(when m ^ n) the rate T^^n^ does not agree with our relation ()28() . Hence, collision rate as 
given by Rautian and Shalagin does not have the property This relation, in turn, is 
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essential in our further discussion. We conclude, that the results of Rautian and Shalagin 
must be treated with care due to quite probable misprints. On the other hand, definition 
presented in this work agrees with expressions used by other authors, see for example |lOl |. 
Eq. 2.4. 
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